Quantum Description of Shuttle Instability 
in Nanoelectromechanical Single Electron Transistor 
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The voltage dependence of nanoelectromechanical effects in a system where the quantized me- 
chanical vibrations of a quantum dot are coupled to coherent tunneling of electrons through a single 
level in the dot is studied. It is found that there are two different regimes depending on the value 
of an applied voltage. If the bias voltage is below a certain threshold value, then the state of the 
mechanical subsystem is located near its ground state. If the bias voltage is above the threshold 
value then the system becomes unstable which manifests itself in the expectation value of the dis- 
placement being an oscillating function of time with an exponentially increasing amplitude. This 
can be interpreted as a shuttle instability in a quantum regime. 
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Nanoelectromechanical systems (NEMS), where electronic and mechanical degrees of freedom are coupled is a 
new and fast growing branch of condensed matter physicsii^. One such a system, a single-electron transistor (SET) 
where the conducting island has a vibrational degree of freedom associated with its center-of-mass motion (Fig.l), 
a so-called nanoelectromechanical-SET (NEM-SET), has been attracting a great deal of attention recently, both 
theoreticallyVv^i^i7A"ii"iiiii2,i.3,i4,i5,i6,i7 experimentallj.l2-La. 
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FIG. 1: Model shuttle system consisting of a movable conducting island placed between two leads, 
acting on the dot from the leads is described by the parabolic potential. 



An effective elastic force 



In Ref. 3 it was shown, that the metallic grain placed between the two leads of NEM-SET becomes unstable and the 
periodic mechanical motion of the grain develops if a large bias voltage is applied between the leads. This phenomenon 
is usually referred to as a shuttle instability (see the review Ref. l20|) . A theory of the shuttle instability developed in 
Ref.'3 was based on the assumptions that both the charge on the island and its trajectory are well defined quantities. 

When we decrease the island size two different quantum effects manifest themself. Firstly, the electron energy level 
spacing in a nanometer-size grain is of the order of 10 K and resonant tunneling effects become essential at small 
enough temperatures. Therefore the single-electron energy spectrum can not be treated as continuous any more as it 
was done in Ref. 0. In this case the characteristic de Broglie wave length associated with the island can still be much 
shorter than the length scale of the spatial variations of the "mechanical" potential. If so, the motion of the island 
can be treated classically. This regime has been studied theoretically in Ref. 0. 

Diminishing the size of the island further (down to 7 A in diameter, for Cgo molecule) results in the quantization 
of the mechanical motion of the island. A NEM-SET system in the regime of quantized mechanical motion of the 
central island was studied theoretically in Ref. It was assumed that the phase breaking processes are strong enough 
to make the density matrix diagonal in the representation of the eigenstates of the quantum oscillator Hamiltonian 
which describes the mechanical subsystem. At the same time it is well known that the expectation value of the 
displacement operator in the eigenstates of a quantum oscillator are zero while it is the coherent state, which is a 
coherent superposition of all |n), that approaches a description of the classical periodic motion of the oscillator as 
h 0. The interesting question arises, therefore, how the expectation value of displacement operator evolves in time 
and what state results when the formal condition of the shuttle instability are satisfied for the quantum NEM-SET. 

The approach used in Ref. |5j was further extended in Ref. to include the position dependence of the tunneling 
matrix elements. However, the possibility of the shuttle instability effects was excluded by assuming a strong external 



2 



dissipation. The case of the weak dissipation was studied in Ref. I6lll2l in the hmit of large voltages. The effect of a 
quantized vibrational mode on electron tunneling through a similar structure (a chain of three quantum dots) has 
been studied in Ref. 0. 

In this article we will show that if a bias voltage between the leads is higher then a certain value, then the quantum 
state of the central island of the NEM-SET evolves in such a way that the expectation value oscillates in time with 
an increasing amplitude. This results shows that the shuttle instability is a fundamental phenomenon which exists 
even when the trajectory of the island and the charge on it are no longer well-defined. 

We use the following Hamiltonian to model our system 

p2 ]\,f 2 v2 

H = Y. ^^>^^L^o.k + e,{X)c^c + — + + T^iXKAic + cU^) . (1) 

a,k a 

Here Aa = J2k ^^k, 

e4X)^eo-dX, (2) 

where d is proportional to an electric field between the leads, 

Ti.flXX) =Toexp{TX/A}, (3) 

where A is a characteristic tunneling length. The first term in the Hamiltonian describes the electrons in the leads, 
the second term relates to the single energy level in the central island, the third and forth terms to the quantized 
vibrational degree of freedom associated with center-of-mass motion of the central island and the last term describes 
tunneling between the electrodes and the island. All energies are measured from the Fermi energy of the leads. We 
assume that the electrons in each electrode are non-interacting with a constant density of states and that all relevant 
energies are small compared to the level spacing in the central island which for typical systems under consideration 
exceeds 100 meV. In this case only one single level in the central island is relevant to the problem. 

It is convenient to introduce dimensionless operators for displacement, x = X/vq, and momentum, p = roP/h, 
where rg = ^h/ (Mwq), and then measure all lengths in units of rp and all energies in units of hwQ. 

The Hamiltonian can now be written as 

H = H,i + Hdot + Ht , (4) 

where 

Hdot = ^[p^ + x^] +€d{x)c^c, (5) 

Hel = ^ eakClak^ak , (6) 

Ht = J2 K^" + ^"^"] ' = ^"(^) ^ ■ 

a 

We will use an equation of motion (EOM) approach to study the time evolution of the operator x in the Heisenberg 
picture. In this approach we deal with the algebra of the time-dependent operators a^kit), c(t), x(t), p{t) in the 
Heisenberg representation whose evolution can be described by a coupled system of EOMs. The EOM for the 
operator aak{t) can be formally solved and eliminated from the system. The reduced system of EOMs is local in 
time under the condition of the wide-band approximation (T>a = const, where Va is a density of state in the lead a). 
It is shown in Appendix 1X1 that the evolution of an operator S = S{x,p,c^c) is described by the following quantum 
Langevin equation 

s = -i[s, Hdot] - J2 " yC-t } [S, Ca] - J2^s, Ci] {zVv^a^ + ^Co,} , (8) 

a a 

where 

a„(i) = Ve-^^°^'a„fc(0), 77, = 2^I?„ . (9) 
^k 

For complete description we need also to know the evolution of c(t), which as shown in Appendix IbI is given by 

cit) = J2 e'''^''B^k{t , 0)a„fe(0) , (10) 

a.k 
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where 



i r* dt 

le -^'i 



(11) 



T is time-ordering operator, f (i) = X^q^o^qI^)' ^a(^) = '7q(3^(^))i ^d{t) = ed{x{t)) and is the operator x in the 
Heisenberg picture. 

We substitute x and p in Eq. © and get 



X + X = F , 



where 



F = dc^c + 2A-1 Y^{~IY ^f^Re [ale] 



(12) 

(13) 
(14) 



To study the shuttle instability in the quantum regime we assume that dimensionless parameters 1/A and d are 
small which allows us to linearize the problem with respect to x. Eq. 1)1 0|l then reads 



a, A: 



dGl 



i-r 



Here = [e.^ - eq + «r/2]-i, E = Vo^T^ To is defined in Eq. @, 

X^k{t) = /*dre*[^°'=-^°+'5l(*-")x(T) 

and the linearized force operator has the following form 

F = Fi + F2 , 

where 

Fi-T^^e'^'""-''' 



'MjO)a^,(0){dG-,G+ +A-1 



(-)"G+ +(-)''G-J} 



is the term of the first order in the small parameters d and A ^ and 



A^ro2^e^('--^')*ai,(0)a^,(0) 



dG- 



i-r 



A 



dG 



(15) 

(16) 

(17) 
(18) 

(19) 



is the second order term in d and A"'^. Thus the above linearization of the force F corresponds to its expansion up 
to second order in d and A~^. 

Averaging Eq. (|13|) with respect to the initial density operator we obtain EOM for the expectation value of the 
displacement of the dot. If we only take into account effects up to second order in Tq (Born approximation), we can 
make the following factorization in the RHS of Eq. 

(20) 



!^2;(*)«L(0)a/3<j(0)) ~ (xit)) fakSak,f3q , 

where Uk = [1 + e(^°'--T'=^/2)/('='^)]-i. Then 

^ + x ^ F [x{t)] , 

where x^t) = {x{t)) ~ x^, xo = {d/2) /a(eo), U{E) = [1 + e^ETeV/2)/ikT)^-i 



F [x{t)] = -2T2 J2 /afelm [Xck] 



dG 



ak 



i-r 

A 



(21) 



(22) 
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It is remarkable that Eq. (|21|l with F given by Eq. (|22(l is the same as the Newton's equation used in Ref. 0. It 
is well-known that for an isolated quantum harmonic oscillator the equation of motion for the expectation values 
of displacement and momentum is exactly the same as the classical Hamiltonian equations (Ehrenfest theorem). 
However, it is not obvious that this statement is valid for the harmonic oscillator coupled to some degrees of freedom. 
Therefore, it follows from Eq. (|22|l that the Ehrenfest theorem is valid for the system under consideration if the 
tunneling is weak. 

We solve Eq. 1)21(1 by Laplace transforms and obtain 



where 



-IE 



xit) 



-XQe^* cost , 



+1 



fa, 



(23) 



(24) 



and fa.±i = /a(eo i 1). The sign of r depends on the bias voltage and can be easily analyzed if the temperature 
is zero. If the bias voltage is below the threshold value Vc, defined by eVc — 2(eo + 1), then fa,+i — and r < 0. 
In this case x{t) exponentially goes to zero (point a; = is stable). If the voltage is above the threshold value, then 
fL,±i = 1, /fl,±i = and r = Td/{2X). In this case x(t) oscillates with exponentially growing amplitude (point x ^ 
is unstable) and the corresponding state of the mechanical subsystem moves further and further from the ground 
state. 

In conclusion, we have studied stability of a system where the quantized mechanical vibrations of a quantum dot 
are coupled to coherent tunneling of electrons through a single level in the dot. Two different regimes have been 
found. For bias voltages above a certain threshold value (which, in the case of zero temperature, is defined by the 
position of the level in the dot plus the energy of a vibrational quantum) the system is unstable. In this regime the 
expectation value of the displacement is an oscillating function of time with an exponentially increasing amplitude, 
which is the signature of a shuttle instability in a quantum regime. Below the threshold the system is stable even 
without any external damping of the dot's motion. 

The author would like to thank L. Y. Gorelik and R. I. Shekhter for valuable discussions. 



APPENDIX A 



The equation of motion for aak(t) is given by 

We formally solve Eq. (|A1|I to obtain 

aofe(t) = *a„fe(0) - i f diie-"-(*-*^)C„(ti) 



Then by using the wide-band approximation (Pq, — const) we get 



where 



If S" = S{x,p,c'c), then its equation of motion is 

a 

Combining Eqs. I|A3|I and (|X5|l yields Eq. 



(Al) 



(A2) 



(A3) 



(A4) 



(A5) 



APPENDIX B 

The equation of motion for c{t) reads 

c(t) = i[H,c] = -tid{t)c{t) ~ iY,Ta{t)Aa,{t) , (Bl) 

a 

where id{t) = ^d{x{t)), Ta{t) = Ta{x{t)) and x{t) is operator x in the Heisenberg picture. Combining Eqs. (jA3(l and 
(|B1|I gives 



c{i) 



edit) -ir^ ' 
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The sohition of Eq. for t > || f H^Ms given by Eq. ^ 



c{t)~iY,y/v^W)a^{t)- (B2) 
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